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Atomic gas far away from the Virgo cluster core galaxy 

NGC 4388 

A possible link to isolated star formation in the Virgo cluster? 



m 
o 
o 

(N 

Ph. 
< 
00 



> 

m 

o 
cn 
o 



B. VoUmer, W. Huchtmeier 

Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 
Received / Accepted 

Abstract. We have discovered 6 lO'^ M0 of atomic gas at a projected distance greater than 4' (20 kpc) from the 
highly inclined Virgo spiral galaxy NGC 4388. This gas is most probably connected to the very extended Ha 
plume detected by Yoshida et al. (2002). Its mass makes a nuclear outflow and its radial velocity a minor merger 
as the origin of the atomic and ionized gas very unlikely. A numerical ram pressure simulation can account for the 
observed Hi spectrum and the morphology of the Ha plume. An additional outflow mechanism is still needed to 
reproduce the velocity field of the inner Ha plume. The extraplanar compact Hll region recently found by Gerhard 
et al. (2002) can be explained as a stripped gas cloud that collapsed and decoupled from the ram pressure wind 
due to its increased surface density. The star-forming cloud is now falling back onto the galaxy. 



Key words. Galaxies: individual: NGC 4388 
dynamics 



Galaxies: interactions - Galaxies: ISM - Galaxies: kinematics and 



Oh 
I 

O 



X 



1. Introduction 

The Virgo cluster spiral galaxy NGC 4388 is located at 
a projected distance of 1.3° from the Virgo cluster cen- 
ter (M87). Its high radial velocity of --1400 kms^^ with 
respect to the cluster mean and its estimated line-of- 
sight distance (TuUy-Fisher method; Yasuda et al. 1997) 
to the cluster center place its three dimensional loca- 
tion very close to M87. In addition, NGC 4388 repre- 
sents one of the nearest Seyfert 2 galaxies and the first 
AGN to be found in the Virgo cluster (Phillips & Malin 
1982). Veilleux et al. (1999) made observations with a 
Fabry-Perot Interferometer in the Ha and 0[iii] A5007 
lines. They found a large, high-ionization plume extend- 
ing northeastwards from the nucleus up to a projected dis- 
tance of 4 kpc^ above the plane of the galaxy. The plume 
is blueshifted with respect to the galaxy's systemic veloc- 
ity. Yoshida et al. (2002) discovered a very large Ha plume 
that extends up to ^^35 kpc northeastwards. This region 
contains ~ 10^ Mq of ionized gas. They argued that the 
dominating source of ionization is the radiation coming 
from the nucleus. Veilleux et al. (1999) discussed 4 differ- 
ent origins for the 4 kpc plume: (i) a minor merger, (ii) a 
galactic wind, (iii) nuclear outflow, and (iv) ram pressure 
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^ We adopt a distance of 17 Mpc to the Virgo cluster. 



stripping. Based on their data they favoured a combina- 
tion of (iii) and (iv). Yoshida et al. (2002) on the other 
hand favoured scenario (i) and (iv). A surprisingly new 
result is that Gerhard et al. (2002) found an isolated com- 
pact Hii region at a projected distance of 17 kpc north and 
4.4 kpc west of the galaxy center. The age of the dozen O 
stars embedded in this region (^3 Myr) implies that they 
have formed far outside the main star-forming regions of 
the galaxy and perhaps even within the hot intracluster 
medium. Gerhard et al. (2002) argued that a radial veloc- 
ity of 150 kms~^ excludes that the involved gas has been 
pushed to this large distance by ram pressure stripping. In 
addition to the studied isolated compact Hii region they 
found 17 Hii region candidates in the data of Yoshida et al. 
(2002). In this article we report the detection of atomic 
gas at a projected distance of more than 20 kpc north- 
eastwards from the center of NGC 4388. The observations 
and their results are described in Sect. Eland Sect. 01 The 
dynamical model is presented in Sect. ^ We compare a 
model snapshot with our observations and discuss possi- 
ble implications in Sect. ]E\ The conclusions are given in 
Sect. El 

2. Observations 

On March, 14-21 2002 we performed 21-cm line observa- 
tions at 5 different positions at a frequency corresponding 



2 VoUmer & Huchtmeier: Atomic gas far away from NGC 4388 

Table 1. Integration times and rms. NGC 4388C 



position 


Center Northwest 


Southwest 


Southeast 


Northeast 


At (min) 


190 80 


120 


120 


350 


rms (mjy) 


1.2 1.7 


1.1 


1.2 


0.9 



to the systemic velocity of NGC 4388 with a bandwidth 
of 12.5 MHz. The two channel receiver had a system noise 
of '^30 K. The 1024 channel autocorrelator was split into 
four parts with 256 channels, yielding a channel separation 
of ~10 kms~^. We further binned the channels to obtain 
a final channel separation of ^^20 kms~^. One central po- 
sition and four positions at a distance of one beamsize 
(9.3') to the NW, SW, SE, and NE from the galaxy cen- 
ter were observed in on-off mode (5 min on source, 5 min 
off source). The integration times are shown in Tabled 
Care was taken to avoid other Virgo galaxies with veloc- 
ities within our bandwidth in all 5 on and off source po- 
sitions. We used 3C286 for pointing and flux calibration. 
The observation times ranged between 80 min and 350 min 
per position. The resulting noise (Table^ is partly deter- 
mined by small amplitude interferences, but is close to the 
theoretical noise of 2 mJy per hour of integration. 



3. Results 

The Effelsberg 100-m spectrum of the central position can 
be seen in Fig. ^ (solid line) . It shows a double- horn struc- 
ture with a '~25% higher peak emission at the receding 
side. The total flux is S'hi=3.9 Jykms"^, which corre- 
sponds to a total Hi mass of A/hi = 2.7 10^ Mq. 

In order to compare our off-center positions to inter- 
ferometric data where the galaxy is spatially resolved, we 
use VLA 21 cm C array data (Cayatte et al. 1990). These 
data have a spatial resolution of 23" x 17" and a channel 
separation of 20 kms^^. We clipped the data cube at a 
level of 4 mJy/beam and cut out the inner 4' x 3' around 
the galaxy center. The resulting spectrum for the central 
position is represented by the dashed line in Fig. ^ The 
total flux and peak flux density of Cayatte et al. (1990) 
are ~20% higher than our values, which is within the cal- 
ibration accuracy. Both spectra agree reasonably well. 

For the 4 off-center positions we have synthesized spec- 
tra from the VLA data using the position of the point- 
ing and the beamsize of the Effelsberg 100-m telescope at 
21 cm (9.3'). We use the spatial information of the VLA 
data to construct the synthesized spectra. Fig.[21shows the 
100-m spectra of these positions (solid lines) together with 
the synthesized VLA spectra. The VLA spectra represent 
Hi gas located in the disk (see Cayatte et al. 1990). If 
there is an emission excess in the off-center 100-m spectra 
this must be due to gas located outside the disk. The rms 
noise of the north-eastern (NE) position is certainly over- 
estimated, because of a relatively strong standing wave 






velocity (km/s) 



Fig. 1. Solid line: Effelsberg 100-m spectrum of the central 
position. Dashed line: spectrum of the VLA data. Dotted 
line: 3(7 noise level of the 100-m spectrum. Heliocentric 
velocities are given relative to the systemic velocity of 
NGC 4388 (wsys=2524 kms^^). 



that had to be subtracted. We made 3 independent ob- 
servations of this position in 3 different nights. In all 3 
sub-spectra the signal is detected. One of them shows a 
flat baseline without a standing wave. The peak emission 
of none of the synthesized VLA spectra exceeds the 100- 
m 3(7 noise level, except the NE position where we have 
overestimated the rms. 

Only the 100-m spectrum of the NE position shows a 
detection. The main peak is around zero radial velocity 
with respect to the galaxy. It has a width of ^80 kms^^. 
In addition, we detect a second smaller maximum at 
160 kms^^, which has the same position and width as 
the corresponding emission of the eastern, receding side of 
NGC 4388. However, it contains ~30% more flux than the 
corresponding maximum of the VLA data. The total flux 
of the emission is 5*111=0.9 Jykms^^, which corresponds 
to an atomic gas mass of A/hi = 6 10^ M©. 

4. The model 

Since the model is described in detail in VoUmcr ct al. 
(2001), we summarize only its main features. The N-body 
code consists of two components, a non-coUisional com- 
ponent that simulates the stellar bulge/disk and the dark 
halo, and a coUisional component that simulates the ISM. 

The 20 000 particles of the coUisional component rep- 
resent gas cloud complexes which are evolving in the grav- 
itational potential of the galaxy. 

The total gas mass is M*°* = 3.810^ Mq, which cor- 
responds to the total neutral gas mass before stripping 
assuming an Hi deficiency of 1.1. To each particle a radius 
is attributed depending on its mass. During the disk evolu- 
tion the particles can have inelastic collisions, the outcome 
of which (coalescence, mass exchange, or fragmentation) 
is simplified following Wiegel (1994). This results in an 
effective gas viscosity in the disk. 
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Fig. 2. Solid lines: Effelsberg 100-m spectra of the four off-center positions. Their locations with respect to the galaxy 
center are marked on top of each panel. Dashed line: synthesized VLA spectra which only show Hi disk emission. 
Dotted line: 3a noise levels of the 100-m spectra. Velocities are given relative to the systemic velocity of NGC 4388. 



As the galaxy moves through the ICM, its clouds are 
accelerated by ram pressure. Within the galaxy's iner- 
tial system its clouds are exposed to a wind coming from 
the opposite direction of the galaxy's motion through the 
ICM. The temporal ram pressure profile has the form of a 
Lorentzian, which is realistic for galaxies on highly eccen- 
tric orbits within the Virgo cluster (VoUmer et al. 2001). 
The effect of ram pressure on the clouds is simulated by an 
additional force on the clouds in the wind direction. Only 
clouds which are not protected by other clouds against the 
wind are affected. 

The non-coUisional component consists of 49 125 par- 
ticles, which simulate the galactic halo, bulge, and disk. 
The characteristics of the different galactic components 
are shown in Table [3 The resulting rotation velocity is 
~200 kms~^ and the rotation curve is flat. 

The particle trajectories are integrated using an adap- 
tive timestep for each particle. This method is described 



Table 2. Total mass, number of particles N, particle mass 
M, and smoothing length I for the different galactic com- 
ponents. 



component 


Mtot (Mq) 


TV 


M (Mo) 


I (pc) 


halo 


2.410'' 


16384 


1.410' 


1200 


bulge 


8.2 lO'' 


16384 


5.0 lO'^ 


180 


disk 


4.1 10'° 


16384 


2.5 10** 


240 



in Springel et al. (2001). The following criterion for an 
individual timestep is applied: 



where is the acceleration of the particle i. The minimum 
of all ti defines the global timestep used for the Burlisch- 
Stoer integrator that integrates the coUisional component. 
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ues with the same gradient, the other part shows a much 
steeper gradient. The minimum velocities are -600 kms~^. 
The minimum absolute velocities are l^l '^800 kms~^. The 
main gas mass is found at a projected distances between 
20 kpc and 30 kpc and radial velocities between -50 km 
and kms^^. At this timestep, the galaxy has an abso- 
lute velocity of ^2000 kms~^ with respect to the cluster 
mean velocity. From Pmax = 5000 cm~^(km/s)^ we derive 
a maximum ICM density of ~ 3 10~^ cm~^. Using 

the ICM density profile of Schindler et al. (1999), this cor- 
responds to an impact parameter of 6 ^ 100 kpc, which 
is consistent with the dynamical model of VoUmer et al. 
(2001). We compare the final timestep {t ~ 120 Myr) with 
our observations. 



Fig. 5. ISM radial velocity as a function of the projected 
distance to the galaxy center. 

The galaxy is on an eccentric orbit within the cluster. 
The temporal ram pressure profile can be described by: 



where inw is the width of the profile (VoUmer et al. 2001). 
We set pmax=5000 cm^3(]^^^/g)2 tHW=50 Myr. The 
efficiency of ram pressure also depends on the inclination 
angle i between the galactic disk and the orbital plane 
(VoUmer et al. 2001). We set i=45°. 

Fig. and Fig. 0] show the evolution of the stellar and 
gas disk during the galaxy's passage through the cluster 
core. The galaxy is seen face-on in Fig. 13 In Fig. 0] the 
position and inclination angle of NGC 4388 {PA = 90°, 
i = 85°) are used. The time of closest approach to the 
cluster center, i.e. maximum ram pressure, is f = Myr. 
Since only the gas is affected by ram pressure, the stellar 
disk does not change during the whole simulation. At t 
—80 Myr the external ~5 kpc ring of neutral gas begins to 
be pushed to the north (Fig.^. At maximum ram pressure 
{t = Myr) this northern extraplanar gas shows 2 distinct 
features: (i) an approximately uniform gas layer in the 
north of the disk and (ii) a horizontal gas ridge located to 
the north east. This is where the accelerated gas, which 
was initially rotating against the direction of the galaxy's 
velocity (at (0 kpc, 20 kpc) in Fig.O and already rotated 
~180°, is now again compressed by ram pressure (at (- 
20 kpc, kpc) in Fig.|2J)- These structures expand in the 
following evolution of the neutral gas. For t > 40 Myr a 
counter-rotating arm is forming (see VoUmer et al. 2001), 
which is responsible for the diagonal gas feature in the 
south-west (Fig. 21). 

In order to show the effect of ram pressure on the ra- 
dial velocities, we show the latter as a function of the pro- 
jected distance from the galaxy center for the last snapshot 
{t ~ 120 Myr) in Fig.EJ From to 30 kpc the radial ve- 
locities with respect to the galaxy center show an almost 
constant gradient from 100 kms~^ to -50 kms~^. At a dis- 
tance of ~30 kpc the velocities bifurcate. One part contin- 



5. Discussion 

Despite the fact that the 100-m spectra do not give precise 
informations about the location of the detected Hi gas, we 
can conclude from our data that the component around 
zero radial velocity lies at a distance of at least half a 
beamsize (4.7' or 23 kpc) away from the galaxy center. 
It contains several 10^ solar masses. If it was closer to 
the galaxy one should see it in the VLA synthesized spec- 
trum if it is high column density gas or alternatively in 
the central 100-m spectrum if it was low column density, 
extended gas. We suggest that it is connected to the very 
extended Ha plume discovered by Yoshida et al. (2002). 

To our knowledge this is the first time that atomic gas 
is observed at these large distances from a galaxy in a 
cluster core region. What is the origin of this extraplanar 
gas? 

Veilleux et al. (1999) discussed 4 different mechanisms 
(cf. Sect. ^ and favoured a mixture of nuclear outflow 
and ram pressure stripping. Yoshida et al. (2002) prefered 
a minor merger and/or ram pressure stripping. If the de- 
tected atomic gas is connected to the very extended Ha 
plume, our data rule out a nuclear outflow as the ori- 
gin of this plume, because of its large total gas mass. A 
superwind that has recently ceased could have provided 
enough energy (see e.g. Heckman et al. 1990), but there is 
no indication for a recent, strong starburst in the center 
of NGC 4388. If the atomic gas and the very extended 
Ha plume represent the debris of a minor merger, its gas 
velocity around relative zero places its orbit in the plane 
of the sky. In our understanding this is inconsistent with 
the morphology of the Ha plume. If the Ha emission rep- 
resents the debris of a tidally disrupted dwarf galaxy, it 
should trace its trajectory. In this picture the curvature of 
the north-eastern part of the very extended Ha emission is 
dynamically difficult to explain. In addition, the Ha plume 
is larger close to the disk, which is also in contradiction 
to the tidal debris scenario. 

We suggest a model, where the atomic gas and the Ha 
are connected and both represent ram pressure stripped 
gas in different phases. In Fig. |S1 (left panel) we show a pos- 
sible ram pressure stripping scenario. The galaxy moves 
to the south and away from the observer. It is on a ra- 
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Fig. 3. Evolution of the model stellar and gas disk. The disk is seen face-on and rotates counter-clockwise. The arrow 
indicates the direction of ram pressure, i.e. it is opposite to the galaxy's velocity vector. The size of the arrow is 
proportional to pv'^^^. The galaxy passes the cluster core at Myr. The time to the core passage is marked at the top 
of the panels. 



dial orbit and its core passage is ~120 Myr ago (Sect.^. 
The galaxy has lost ~90% of its atomic gas. We show 
only gas with a local density greater than n ~ 30 cm~^. 
For clarity the greyscales are inverted, i.e. dark regions 
correspond to low column density gas. Since the stars 
are not affected by ram pressure, the stellar disk is sym- 
metric. The gas disk is strongly truncated at a radius of 
^5 kpc (~ 1'). This is in agreement with the VLA ob- 
servations of Cayatte et al. (1990). Only the high density 
eastern edge of the stripped, extraplanar gas distribution 



can be seen. In addition, the counter-rotating, backfalling 
arm in the south-west has high local densities and thus 
appears in Fig. El (left panel). The main off-plane gas 
mass is located at (-20 kpc, 20 kpc). For simulations with 
smaller maximum ram pressure, it is located closer to the 
galaxy. The corresponding velocity field is shown in Fig. El 
(right panel). Within the disk plane, the velocity field is 
dominated by rotation. The western side of the disk is 
approaching the observer. The stripped gas at distances 
smaller than 10 kpc above the disk plane has positive ve- 
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Fig. 4. Evolution of the model stellar and gas disk. The position and inclination angle of NGC 4388 {PA — 90°, 
i = 85°) are used. The arrow indicates the direction of ram pressure, i.e. it is opposite to the galaxy's velocity vector. 
The size of the arrow is proportional to pfg^i- The galaxy passes the cluster core at Myr. The time to the core 
passage is marked at the top of the panels. 



locities, whereas the gas at distances greater than 10 kpc 
has negative velocities. The counter-rotating, backfalling 
arm has a mean radial velocity of ~ kms^^. 

Fig. [7| shows the model 100-m spectrum at central 
position (left panel) and the northeastern (NE) position 
(right panel). We find a total Hi mass of ~ 5 10^ M© and 
~ 7 10^ Mq at the central and north-eastern position, re- 
spectively. 

The central spectrum shows a double horn structure 
and the receding side has a higher peak flux as it is ob- 



served (Fig.^. However, a part of the flux from the high 
velocity end of the approaching side is missing, mainly 
because it has a lower local density. 

For the north-eastern position we observe a clear re- 
semblance between the model and the observed spectrum 
(Fig. 121). The maximum at 160 kms~^ is due to the ex- 
tended eastern, receding side, whereas the main maximum 
around '-~^-50 kms^^ radial velocity is due to the mate- 
rial at distances of about 20 kpc from the galactic disk 
(Fig. EI) . Thus the ram pressure scenario can account for 
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Fig. 6. Left panel: stellar distribution as contours, gas distribution as greyscales. Only gas with a local density greater 
than ~ 30 cm~^ is shown. For clarity the grayscales are inverted, i.e. dark regions correspond to low column density 
gas. Right panel: model velocity field. Dashed lines mark positive, solid lines negative velocity contours. The velocities 
range between -200 kms~^ and 260 kms~^. The contour stepsize is 20 kms~^. 




Fig. 7. Left panel: model 100-m spectrum at the central (C) position. Right panel: model 100-m spectrum at the 
northeastern (NE) position. Velocities are given relative to the systemic velocity of NGC 4388. 



the very extended tail. However, it is not possible to re- 
produce the velocity field of the NE Ha plume (Veilleux 
ct al. 1999), since it shows a blueshift with respect to the 
galaxy's systemic velocity. In order to explain this feature 
an additional outflow mechanism is needed. 

How can this gas survive within such a hostile envi- 
ronment like the intracluster medium (ICM)? Our model 
is not able to investigate this subject, because it considers 
only the neutral gas phase and does not include ISM/ICM 
shocks that might heat the gas. Evaporation and ioniza- 
tion physics are also not included. Thus we have to rely 
on analytical arguments. 

Veilleux et al. (1999) argued that in the case classical 
thermal heat conduction a typical Hi cloud evaporates in 
the hot ICM within 10^ yr. On the other hand, Vollmer 



et al. (2001) showed that within the ICM evaporation is 
saturated. If the stripped atomic gas clouds have column 
densities of 10^"'^ cm~^, the evaporation time can reach 
~10* yr. They also discussed the possibility that the dense 
part of the stripped neutral gas contracts due to its self- 
gravity, because the internal heating (stellar radiation and 
supernova shocks) are absent. Since the dust is stripped 
together with the gas, molecules might form when the gas 
density becomes high enough. The only gas heating is then 
the X-ray radiation from the hot ICM. An equilibrium 
cloud has a mass of several thousand solar masses, a ra- 
dius of ^10 pc, and a column density of a few 10^^ cm~^. 
Only 10% of this column density is in form of neutral 
hydrogen. Could this be the case for NGC 4388? These 
clouds can in principle survive long enough to be pushed 
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to the observed large distances. It is then not excluded 
that these clouds collapse and form stars. 

Interestingly, Gerhard ct al. (2002) found a compact 
Hll region that has formed stars in situ, far away from 
the main star-forming disk. This Hll region is located at 
a projected distance of 17.5 Mpc from the galaxy center 
in the northwest of the galactic center and has a radial 
velocity of 150 kms~^. Within our model it is possible 
that this Hii region is falling back to the galaxy after 
having been pushed by ram pressure to a large galacto- 
centric distance. When a stripped cloud collapses, its col- 
umn density increases rapidly and it decouples from the 
wind (ram pressure). It then begins to fall back to the 
galaxy due to the galaxy's gravitational potential. With 
the relatively low metallicity (Gerhard et al. 2002) this 
stripped gas was located approximately at the optical ra- 
dius (i?25) of NGC 4388 before stripping has taken place 
(Due et al. 2000). This scenario is in agreement with the 
positive radial velocity of the compact Hii region. In ad- 
dition, Gerhard et al. (2002) found a number of Hll region 
candidates in the very extended Ha plume. Within our 
picture this would mean that the stripped gas clouds have 
become molecular and are now forming stars. This mecha- 
nism woiild be analog to the formation of globular clusters 
in the tidal tails of gravitationally perturbed galaxies. 



to thank V. Cayatte and J. van Gorkom for providing us 
kindly their VLA data. 
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6. Conclusions 

We have observed the Virgo cluster core galaxy NGC 4388 
in the 21 cm Hi line with the Effelsberg 100-m telescope 
at 5 different positions: one in the center and 4 located 
one 100-m beamsize to the NE, NW, SW, and SE from 
the galaxy center. The results of these observations are 
compared to a realistic dynamical model including ram 
pressure stripping. The results are: 

1. We detected extraplanar Hi gas of ^6 10'' Mq located 
at a distance of at least 20 kpc to the north-east of the 

galaxy center. This gas is most probably connected to 
the very extended Ha plume detected by Yoshida et 
al. (2002). 

2. The dynamical model is able to reproduce the Hi defi- 
ciency, the truncated Hi disk and the extraplanar emis- 
sion. NGC 4388 passed the cluster core ~120 Myr ago. 

3. Within this scenario the maximum ram pressure is 
Pmax = 5000 cm~'^(km/s)~^, the maximum ICM 
density encountered during core passage is nicM ~ 
3 10~'^ cm^''; the galaxy's absolute velocity with re- 
spect to the cluster mean is ~2000 kms~^, and the 
angle between the disk and the orbital plane isi = 45° 

4. The extraplanar compact Hll region recently found by 
Gerhard et al. (2002) can be explained as a stripped 
gas cloud that collapsed and decoupled from the ram 
pressure wind due to its increased surface density. 
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